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Polar Factorization of a Matrix"" 

John Z. Hearon** 

(April 24, 1967) 

It is known that if ^ is a bounded linear operator with closed range on a Hilbert space then A can 
be factored as A = UH, with U a partial isometry and H nonnegative and self adjoint. For the finite- 
dimensional case a strictly matrix-theoretic derivation is given based on the concept ol a general- 
ized inverse. Certain properties of the factors are given as well as conditions under which H or both 
U and H are uniquely determined by A. A pivotal item in the derivation is the representation of a square 
partial isometry as the product of a unitary matrix and an orthogonal projection. This representation 
is new, of some interest in itself and greatly simplifies the derivations. 
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1. Introduction 

It is well known but not well discussed in the 
matrix literature' that a square matrix, A, can be fac- 
tored as ^ = UH where C/ is a partial isometry and 
H is positive semidefinite. The primary purpose of 
this paper is to give a fairly direct demonstration of 
this factorization, and the main result is thus not 
new (see footnote 1). The demonstration, as well as 
the deduction of certain properties of the factors, is 
based on a characterization of partial isometrics which 
is new and of some interest per se. While not strictly 
necessary it is possible and illuminating to cast part 
of the development in terms of generalized inverses 
of singular matrices. 



2. Notation and Preliminaries 

In what follows all matrices are considered to have 
complex entries. We denote by p{A), R(A), N{A) and 
y4* rank, range, null space and conjugate transpose, 
respectively, of any given matrix. When A is non- 
singular, A~'^ denotes the inverse. For generahzed 
inverses a special terminology is used. This termi- 
nology, previously introduced and related to others 
[5, 6]- is as follows: For a given matrix A denote by 
C\{A) the set of all matrices B such that ABA=A. 
Then Co{A) is defined as the set of all matrices B such 
that BeCiiA) and AeCx[B)\ C-iA) is the set of all mat- 
rices B such that BeCiiA) and AB is hermitian; finally 
C4{A) is the set of all matrices B such that BeCs{A) 



*An invited paper. 

**Present address: Mathematical Research Branch. NIAMD. National Institutes of 
health, Bethesda, Md. 20014. 

' The factorization, with condition for hoth factors to be unique, cf. Theorem 2 to follow, 
is tiiven as a problem in [4, p. 171]. The factorization is well known as a result for bounded 
operators with closed range on a Hilbert space [1, 8]. Desoer and Whalen [1] give the fac- 
torization where U has the property that U* is the pseudo-inverse (which on a finite-dimen- 
sional inner-product space is the Moore-Penrose inverse) of U. This is equivalent to f/ being 
a partial isometry. 

2 Figures in brackets indicate the literature references at the end of this paper. 



and BA is hermitian. We note that the set C4{A) con- 
tains a single uniquely determined matrix which is 
the Moore-Penrose generalized inverse [7]. We call 
a matrix BeCiiA) a C /-inverse of A. The relation be- 
tween a C.rinverse, as here defined, and the "weak 
generahzed inverse" of Goldman and Zelen [3] has 
been noted elsewhere [5J. Repeated use will be made 
of the following fact: If BeCiiA) then p{B) ^ p{A) 
= p{AB) = p{BA), with strict equahty if and only if 
BeCiiA) [5, 9J. 

We call a matrix A a partial isometry if there exists 
a subspace, S, such that x'^^A'^Ax = x'^x, when xeS, 
and Ax = 0, when xeS^, where S^ is the orthogonal 
complement of S. This dehnition is equivalent to the 
requirement that A'^A be an orthogonal projection [2], 
[4, p. 150J. 

3. The Polar Factorization 

We begin with the following two lemmas 

Lemma 1. The square matrix, A, is a partial isom- 
etry if and only i/A = QE where Q is an isometry and 
E is an orthogonal projection. 

Proof. If A = QE, with (?*(?-/ and E = E^ = E'', 
we have A'^A = E and A is a partial isometry. Let A 
= QH be the usual polar factorization of A, where Q 
is unitary and //is positive semidefinite. If ^ is a partial 
isometry then A^A = H'^ is hermitian and idempotent. 
If so then //, the positive semidefinite square root of 
//^, is also hermitian idempotent. 

Remark. It is an obvious consequence of Lemma 1 
that ^ is a partial isometry if and only if A = FQ where 
F is an orthogonal projection and Q is unitary. For, 
from A = QE we have A = QEQ'^Q, and we identify 
F with the orthogonal projection QEQ^. Conversely 
A=FQ = QQ''FQ = QE. 

Lemma 2. Let A be normal and BeCi{A). Then if 
E = AB is normal, E is uniquely determined by A, 
and EA =AE = A, 
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Proof. From ABA=A = EA it follows that Ax = kx, 
implies Ex = x, provided X 7^ 0. Let p{A)=r. Then 
there are linearly independent xi such that Exi = Xi, 
1 ^ i ^ r, and since p(E)=p {A) we have R{E)=R {A) . 
Since £" is a normal projection, it is an orthogonal pro- 
jection and thus uniquely determined by its range and 
hence by A. Further, since E and A are normal, E"^ = A 
shows that N{E) =N{A). Hence E and A have a com- 
plete set of eigenvectors in common and must commute. 

Theorem 1. Let A be any n-square matrix. Then 
there exists a partial isometry U and a positive semi- 
definite matrix H, such that 

(i) A = UH, 

(ii) U*A = H, 

(iii) N(U) = N(H)-]NiA), 

(iv) U maps all of w- space onto R(A) 

(v) UH = HIJ if and only if A is normal, and in this 
case U is normoJ. 

Proof. Let A = QH be the usual polar factorization 
of y4, where Q is unitary and H is positive semidefi- 
nite. Let P be any Ci-inverse of H such that E = HP 
is an orthogonal projection. Then HPH = EH = H 
and we hdve A = QH=QEH= UH. where U=QE is, 
by Lemma 1, a partial isometry. Thus (i) is proved. 
Now U''U = E and hence, from (i), U''A=EH = H 
which is (ii). It is clear that N(U) = N{E) and that 
N{A) = N{H): and since p{E)=-p{H), £'// = /7 shows that 
N{E) = N{H). Thus N(U) = N(E) = N(H) = N{A), which 
gives (iii). (^iven A, the projection E = HP is, by Lemma 
2, uniauely determined. If P is chosen to be non- 
singular, as is plainly possible (see after (2) below), 
then A = UH=UEP-' = UP-\ and (iv) is evident. 
Suppose A to be normal. Then A = QH = HQ and 
from this and EH = HE = H, which we have from 
Lemma 2 (but which in this case is obvious from 
EH = H since E and H are hermitian), it follows that 
EA = AE = A. Bu\ then AE=A = HQE = HU=UH. Con- 
versely, suppose UH = HU, We have at once that 
AE = EA = A, which shows ihai IV{A^) = N{E) = N (A) . 
Given this, and N{H)=N{E),we have from ^* = //(?* 
that QyeN{E) whenever yeN(E). We can now assert 
that HQy=QHy = 0, when yeN{E). Finally, HU=HQE 
= UH=QH implies that HQx = QHx, when xeR{E). 
We have proved that HQ = QH and hence that A is 
normal. Given this, horn A^QH = QEH = HQ = EHQ 
= EQH, we have QEx = EQx when xeRiH), and we 
have seen that ()yeM£) = M//) when yeN{E). Thus 
QE = EQ and U is normal. 

Theorem 2. Let A = UH, where U is a partial isom- 
etry and H is positive semidefinite. Consider the con- 
ditions: (i) U*A = H, (ii) p(U) = p(H), (iii) N(U)-N(H). 
Then, if{i) holds, H is uniquely determined; (Hi) holds 
if and only if(i) and (ii) hold, and in that case both U 
and H are uniquely determined. 

Proof. By Lemma 1, we may replace Uhy ^E" with 
Q unitary and E an orthogonal projection. Then, if (i) 
holds, U^A = EH = H. This being so, we hdve A = UH 
= QEH=QH, and H^ = A''A. Thus H is the unique 
positive semidefinite square root of A'^A. We next 
show that (iii) is equivalent to (i) and (ii) together. 
Let (i) and (ii) hold. Then, with U=QE, (i) gives C/*^ 
= EH=H, which with (ii) impHes (iii). Let (iii) hold. 



We obviously then have (ii). Further, with U=QE, 
(iii) states N{E) =N{H). Let xi, X2, . . ., Xr be any 
orthonormal basis of N{E) =N(H). Then from E = I 
— Xxixf, we have EH = HE = H. This being the case, 
U''A=EH = H which is (i). Now let A = UxHx = UiH2 
be any two factorizations of A and assume (iii). Since 
(iii) imphes (i), Hi = H2 = H and we have U\H=U2H 
which implies UiX^Uox, when xeR{H). But (iii) now 
also requires N(Ux) = N{U^y) = N{H) and hence Uxy 
= Uoy for yeN{H). Thus U, = U.. 

If H is hermitian, then H=T diag (A, 0)7*, where 
T is unitary, A is real, diagonal and nonsingular. In 
the following discussion let this unitary similarity via 
T be denoted by // — diag (A, 0). Then for arbitrary 
K, L, and D of appropriate sizes and shapes any P 
such that 



A 1 K 

L D 

is a Ci-inverse of H. For, from 

7 AK 



(1) < 



E = HP 







(2) 



we have that HP is idempotent and has the rank of H 
and this is known [5] to be necessary and sufficient 
for PeCi{H). Now E, in (2), is hermitian if and only if 
K = 0. Thus given K = 0, any P as in (1) will serve in 
the proof of (i), (ii), (iii) and (v) of Theorem 1, and any 
P as in (1) with D nonsingular will serve in the proof 
of (iv) of Theorem 1. Now we could, in the proof of 
Theorem 1 except for (iv), forthwith have taken PeC^iH) 
or PeCsiH), for in both cases E = HP is hermitian. For 
the proof of (iv), we could have then noted that for 
PeC^iH), P-\-E{) is nonsingular when Eo is the princi- 
pal idempotent matrix of// (and of P) associated with 
the zero root,"^ and H{P^ Eo) = HP = E. Of course 
the Theorem 1 could be proved, without reference to 
generalized inverses, by simply producing P as in 
in (1) with K=0, noting that E as in (2) is then hermi- 
tian idempotent, and tliat, subject to K = 0, E is in- 
variant under choices of P. The pivotal idea of the 
proof is the observation that given A = QH, we have 
(i) of Theorem 1 at once, in view of Lemma 1, if we 
can produce an orthogonal projection, E, such that 
EH = H. This possibility is suggested by considering 
generalized inverses and that it is indeed possible is 
perceived at once by considering the Moore-Penrose 
generalized inverse, but as we have seen, other "in- 
verses" will serve as well. 

In the proof of (iv) of Theorem 1 and in the above 
discussion we have encountered an observation which 
may be set out as a corollary. 

Corollary. // A is any square matrix, there exist 
matrices P such that AP is a partial isometry. Further 
there exist such matrices P which are normal, in par- 
ticular positive definite. 



the ort.hogonai |)roi«jction upon N{H) =N{P). 
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Proof. As we have seen any P as in (1) with /C=0 
has the required property. Any P as in (1) with K = 0, 
L = 0, D normal and nonsingular is normal and non- 
singular and has the required property. In particular, 
\i K = 0, L = 0, and D is positive definite, we have a 
positive definite P from (1). 

From Theorem 1, the corollary and the usual polar 
factorization A = QH, we have the following statement: 
If A is any square matrix, there exists an isometry Q 
and a partial isometry U such that Q^A = IJ'^A = H, 
where H is positive semidefinite. If A is nonsingular 
there exists a positive definite matrix, C, such that 
AC = Q is an isometry, but there always exists a posi- 
tive definite P such that AP= U is a. partial isometry. 
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